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■ rapid method of analysis for CORRELTITIOH -DT- 
ASSET RETURN TO FUTURE FINANCIAL LI AB^^^^^^^ 

The present invention is related genera U3«ic<wV ^ffiSCfKM • 
and system for selecting a portfolio of assets for achieving 
optimum correlation of asset return to a selected standard 
financial index. More particularly, the invention is related 
to a highly efficient, rapid method and system for choosing an 
asset portfolio having the optimum correlation of the asset 
return to a time dependent financial index, such as a financial 
liability, at each of a number of selectable asset- return 
levels. 

Management o£ portfolios of assets has historically 
emphasized maxlml.m, the return on assets with the objective 
of at least ootperformin, the market. However, in some 
financial industries the objective, or the £i,ure of merit, is 
also related to meeting future liabilities rather than Just 
•chievin, hi,h return on assets. Fre,„ently, an institution 
"Ul have a future intended use of the assets which requires 
their availability at some future time, if assets ere 
Performin, differently than llebillty requirements, substantial 
financial hardship can result. For e„mple, insurance 
companies and corporate pension plans have well defined future 
financial liabilities which must be met. Consequently, 
althou,h return on assets is one important objective, meeting 
future liabilities is also important and can be even more 
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important in many instances. In fact, many pension plan 
managers are now required to meet the standards set forth in 
FASB Statement 87 (Financial Accounting Standards Board) on 
pension fund accounting. Under the FASB Statement a market 
interest rate return on pension funds is the standard index and 
is to be based on A-rated ten year corporate bonds. Under this 
FASB Statement any deficit in corporate pension funds are now 
reflected on the balance sheet. Any such deficit would 
therefore have substantial adverse effect on the apparent net 
worth of the subject corporation. Consequently, this FASB 
Statement standard strongly encourages maintenance of a surplus 
for a pension fund. As an example of the importance of 
matching the liability requirements under the FASB Statement, 
consider the percentage change possible for pension plan 
liabilities, as measured by the accumulated benefit obligation 
(ABO). If, for example, interest rates increase by 1% in one 
year over the present rates, the present value of the ABO would 
decline by 10% if the ABO has a duration of 10 years. 
Likewise, if interest rates were to drop by 1% in one year, the 
present value of the ABO would increase by 10%. The potential 
for such dramatic fluctuations in liabilities clearly deserves 
careful attention by parties obligated to meet future liability 
streams. 

In order to timely meet future financial liabilities 
and maintain a proper surplus fund for a pension plan, a number 
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of methodologies have arisen including "immunization**, "cash 
matching"; and some preliminary efforts have even been directed 
to utilizing stock funds. 

The "immunization" method of meeting future financial 
liabilities uses bonds having substantially the same duration 
as the liability stream. Duration is a measure of volatility 
expressed in years, which is similar to, but more precise than, 
average life. The duration is calculated as the weighted 
average amount of time to the receipt of the payout. There are 
however significant drawbacks to "immunization", with one 
primary disadvantage being the relatively low excess return on 
assets generally achieved by the method. Additional 
limitations are imposed by the two major assumptions made by 
the strategy: The yield curve (a plot of yield to maturity on 
bonds versus their time to maturity) will only make parallel 
shifts. Consequently, regardless of maturity, when market 
conditions change, all bonds allegedly move exactly the same 
amount in yield. This clearly is not the case since there have 
been substantial inconsistencies in the past for the difference 
in rates for short-term bonds and long-term bonds. Secondly, 
all cash flows in excess of required annual payments can 
allegedly be reinvested at the yield to maturity of the 
portfolio. This presumption is also clearly not true since 
sharply declining or rising interest rate environments will 
make it extremely difficult to carry out reinvestment. 
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Furthermore, this strategy does require more ongoing management 
of the portfolio in order to sell or buy more securities to 
match the actuarial schedule and maintain a proper 
asset/liability match. 

The "cash matching" method utilizes a bond portfolio 
having numerous component bonds with various maturity dates and 
payout rates to precisely match the liability requirements of 
the pension plan. Such an approach has the same primary 
disadvantage as the "immunization" method and further requires 
additional effort to assemble the portfolio. Frequently, the 
"cash matching" method demands payment of a premium to achieve 
the correct, mix of bond rates and maturity. Both of the first 
strategies ("immunization" and "cash matching") must invest in 
fixed income securities to provide the assurance of receiving 
the necessary cash" flows. In fact, they must primarily invest 
in U.S. treasury obligations since investments in corporate or 
mortgage securities increase the chance for default or for call 
risks which can have the effect of changing the projected cash 
flow. 

Pension plan liabilities or other future liabilities, 
such as are present in the insurance industry, are long term in 
nature. Therefore, a future liability stream can greatly 
benefit from the compounding effect of investment in higher 
returning assets, such as common stocks. However, attempts to 
characterize stocks in terms of a time duration parameter or 
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otherwise have not been successful. In the last few years many 
unsuccessful attempts have been made to develop a system 
whereby a portfolio of equities is linearly optimized relative 
to a liability stream. There have been attempts to parallel 
the "cash matching" techniques with the use of stocks, instead 
of bonds. This approach has involved matching the expected 
dividend flow of the portfolio to the liability stream. 
Unfortunately, stock dividend yields are unpredictable, 
particularly beyond 3 years in the future. Another major 
effort in equities has been directed to an "immunization" type 
treatment. In this effort an attempt was made to calculate the 
duration of stocks on an individual basis, as well as on a 
portfolio basis; but these attempts also have been 
unsuccessful, primarily due to the long term unpredictability 
of stock dividends. 

In a related patent application, incorporated by 
reference herein and having serial number 281,560 and filed 
December 8, 1988, an improved method and system were set forth 
directed to correlating return on assets to a financial 
objective over time. In performing the analysis to determine 
the optimum assets of a portfolio to track the financial 
objective, the machine time and efficiency of the evaluation 
process can limit the number of assets considered in 
constructing the portfolio. Such limitations on the number of 
assets which are considered for inclusion in the optimum 
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portfolio can also limit the performance of the selected 
portfolio* 

It is therefore an object of the invention to provide 
an improved method and system for determining the optimum 
portfolio of assets for tracking a financial index. 

It is another object of the invention to provide a new 
method and system of efficiently selecting the optimum 
portfolio of assets for tracking a financial index. 

It is ah additional object of the invention to provide 
an improved method and system of rapidly analyzing a large 
number of potential assets to select the optimum portfolio of 
assets to track a, financial index. 

It is a further object to provide a new method and 
system of enlarging the number of potential assets under 
consideration for incluision in a portfolio of assets, while 
reducing the time required to select the portfolio of assets 
which best track the behavior of a financial index. 

It is..,aneither object of the invention to provide a 
rapid/ more efficient method and system of selecting the 
weighted values for assets selected from a universe of possible 
assets for a portf<^io designed to track a financial index. 

It is an -additional object of the invention to provide 
an improved method and system for reinvesting cash flow from a 
portfolio starting with that current portfolio of assets. 
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Further objects and advantages of the present 
invention, together with the organization and manner of 
operation thereof, will become apparent from the following 
description when taken in conjunction with the accompanying 
drawings described hereinbelow. 



BRIEF DESCRIPTION OF THE DRAWINGS 

FIGURE lA is a functional flow chart illustrating 
operation of one method of portfolio construction and FIG. IB 
is a flow chart illustrating operation of a particular method 
of the invention; 

FIGURE 2A illustrates simulation results of annual 
total liability returns of a pension plan (dashed) and annual 
total asset returns (solid line) for a preferred form of the 
invention and FIG, 2B shows annual total liability returns of 
the simulated pension plan (dashed) and annual. total asset 
returns (solid line) for the Standard & Poors 500; 

FIGURE 3A is a bar graph of simulation results for 
funded status returns over time for a preferred method of the 
invention and FIGURE 3B is a bar graph of funded status returns 
over time for the Standard & Poors 500; 

FIGURE 4 a comparative plot of cumulative funding 
status for simulation results over time for a pension plan 
liabilities (dashed), a portfolio derived by a preferred method 
of the invention (solid line) and the Standard & Poors 500 
index (dashed and dotted); 
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FIGURE 5A shows the correlation between asset return 
and liability return for a preferred method of the invention 
and FIG. 5B illustrates correlation between asset return and 
liability return for the Standard & Poors 500 index; and 

FIGURE 6 illustrates the boundary line of minimum risk 
for various future asset return levels. 



DETAILEP DESCRIPTIPW OF PREFERREP EMPOPIMEHTS 
Broadly Cstated, a method and system are described for 
selecting a portfolio of assets and correlating a future asset 
return of the portfolio to a financial index, such as, a 
liability inde^/ an inflation index, or any other accepted 
index and mixtures thereof . Specific examples of indices are 
liability indices, such as, individual corporate pension plan 
liabilities and insurance company liabilities. The consumer 
price index and wage growth index are examples of an inflation 
index, and other indices can include accepted stock price 
indices and futures markets indices. The method includes 
selecting asset portfolios which optimally correlate portfolio 
returns to the future desired payouts or payments needed over 
time to fulfill the desired financial objective. In the 
general case the user selects a standard index to which optimum 
correlation is desired for the selected portfolio having a 
future asset return over time. The process of selecting the 
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standard index can involve obtaining input (such as actuarial) 
in terms of the characteristics of future cash payments 
discounted to present value based on a range of discount rate 
and wage (inflation) values. This information can be used to 
construct a functional behavior for the present value of the 
liability. A decision is then made, such as by a company 
pension fund manager, that certain discount rates and inflation 
assumptions should be made. On this basis the current 
liabilities are projected back in time using these assumptions 
and a plurality of assets are examined to determine their 
sensitivity to the past behavior of the liability returns. In 
the most general sense if one can determine an index to which a 
portfolio of assets has a strong correlation, this sensitivity 
can be used to select a set of assets which will match the 
behavior of the index as it changes over time. As a particular 
example an actuary can provide specific ranges of present value 
liability for a range of discount rates and inflation rates. 
The change over time of the liability from month to month over 
a twenty four month period can yield a liability return. The 
analysis to be described in more detail hereinafter determines 
which selected ones of a plurality of assets track the 
liability returns with best correlation. The resulting 
weighted set of assets form the portfolio to follow the future 
liability returns. An analysis using the selected standard 
index can be performed on a plurality of assets, such as, for 
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example, at least one of the following categories of assets; 
stock securities, real estate investments, futures contracts, 
options, commodities, currencies and precious metals. The 
analysis allows the identification of the combination of weight 
percentages of selected ones of the plurality of assets 
yielding the optimum correlation of the future asset return to 
the standard index. Optimum correlation is thus achieved by 
calculating a minimum standard deviation or a variance for the 
difference between the return of the portfolio of assets and 
the selected standard index return. This method and system are 
particularly applicable for, but not limited to, the insurance 
industry and management of pension fund liabilities* 

FIGS. lA and IB illustrate in functional flow charts 
the procedures followed in carrying out two forms of the 
invention. In the first method shown in FIG. lA (and described 
previously in pending application having serial number 281,560) 
the correlation of the expected asset return of a portfolio to 
a standard index one is initiated by input of various basic 
information. This information includes, for example, 
establishing the fundamental statistical characteristics of 
liability returns, and future payment schedules for matching a 
desired index, such as the future stream of financial 
liabilities of a pension plan. As described hereinbefore, the 
future payment schedule for a pension plan can be determined by 
using actuarial data. These future liabilities can be 
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characterized in terms of an accumulated benefit obligation 
(ABO), that is, the price you would have to pay if the 
liabilities were sold at a selected time. The total outlay 
required to pay retirement wages for the pension plan are 
discounted back to the present value at the market rate 
interest (currently 10%) . other related characterizations can 
be used, such as a projected benefit obligation (PBO), by 
accounting for inflation in the growth of wages at retirement. 
This amount is converted to a percentage and an expected salary 
at retirement, discounted to present value. Therefore, 
although the ABO is affected primarily by interest rates, 
additional standard measures, such as the PBO, account for 
inflation. Therefore, the method is also generally effective 
for calculating the convolution of complex effects with one 
another. The method only requires optimizing correlation of 
the time behavioral performance of future asset return relative 
to the particular standard index, which includes any 
conceivable selected characteristic which assets are found to 
be sensitive to. 

In the manner illustrated in step 1 of FIG. lA, 
various input files are therefore created to begin the 
analysis. These input files can include, for example, asset 
return information for the universe or plurality of assets to 
be sampled in the analysis. Also established as data files are 
the data representative of the standard asset return over time, 
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such as target returns for a future liability stream of a 
pension plan or an insurance company. The future liability 
stream can depend on interest rates and/or inflation rates and 
other variables which can affect the liability stream. For 
example, ds described hereinbefore, a surface can be generated 
which describes the behavior of liability return as a function 
of both intedest rates and inflation rates. Other information 
in the data files can be identification information for the 
plurality of assets, current price and market capitalization of 
the assets, as well as the characteristic weight percentages of 
assets in a previously selected portfolio. Weight percentages, 
for example, from a prior period would be used in the most 
preferred embodiitfent . 

As illustrated in the second box of FIG, lA, the 
number and type of the plurality of assets should be specified, 
the time period for matching the standard, the number of 
returns used to calculate a statistical correlation and the 
number of asset returns used to calculate averages for the 
plurality of securities to be analyzed to select the optimum 
set of weighted assets. 

The method in FIG. lA then initializes information 
preparatory to analyzing the plurality of assets, such as, 
establishing names of securities, associated identifier 
information, industry codes, prices of securities, market 
capitalization, weight and percent of the previously calculated 
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prior portfolio, the period for the asset return, the name of 
the target or standard for measuring a standard asset return 
and the time period for the standard returns. 

In the next box in FIG. lA the correlation between the 
future asset return and the standard index is optimized by- 
first generating a covariance array. While other nonlinear 
statistical analyses are possible, this method being described 
is a preferred method of carrying out the analysis. For 
example, another useful statistical method of analysis is 
correlation parameterization which is embodied in the computer ' 
software program Appendix III. As shown in the next step in 
FIG. lA, the average return of each security is calculated 
followed by imposing certain constraints on the calculation, 
such as setting a range of weight percentages to be tried. The 
calculation is then implemented to a solution by a standard 
computer program quadratic technique (see Appendix I). This 
step is then followed by determination of various statistical 
parameters, such as X and B, standard error, portfolio returns 
over various time periods and for selected weights. The 
analysis is then completed by printing output (see attachment 
to Appendix I) such as asset weights, sensitivity factors for 
selected assets of the portfolio, statistical parameters, 
sorted buy and sell orders and sector weights. 

A simple example of utilizing the preferred 
statistical method is illustrated for a portfolio containing 
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three stocks (designated 1, 1, and 3). In order to find the 
optimum weight percent for each of the three stocks in the 
portfolio, the minimum standard deviation (square root of 
variance) is calculated for the differences between the assets 
of the portfolio and the future liabilities as represented by 
the standard asset return over time. The risk is therefore 
defined as the standard deviation of differences: 

- ^( R . i) - n R R^]^'^^ "^t^n " " variance 



Risk 



1/2 



c rv - I ) 



where: ?pi * total return on the portfolio during period i; 
R.pj « total return on target or standard in period i, 
(T" = standard deviation; 

Rp « average return on portfolio, i » l,n; and 

R^ - average return on target or standard portfolio of 

assets, 

1 * 1, • • • , n*,» • 

The portfolio return equals percentage weight for each stock 

times the return on that stock: 

^ 



R . » 5 x.R . 



1" 

B the weight in the portfolio of the stock j 

Rgj J - the return on stock j in periods i 

Now in the definition oiE risk, as set forth above, we can 

substitute the following: 
3 
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Making this substitution, a determination of risk in the manner 
set forth above results in the calculation of the covariance of 
the stock with each of the other stocks in the portfolio after 
subtracting the return of the target, or standard index, from 
the future asset return of each stock. 

The covariance of stock 1 with stock 2 is therefore: 

We calculate all the spanning covariances and put them 
in a matrix form: 



COV [(B^ - R^), (Rj - R^)] COV [(Rj - R^), (R^ . R^)]" 

COV [(Rj . R^), (R^ - R^)) Cov [(R3 - Rj), (Rg . R^)] 



In order to calculate the risk, we add up all the Cov 
terms times the weights in each stock: 



Risk = r(Rp . R^, Jt Z Z h^, cov ((R^ - R^), (R . 8 )jj 



In order to minimize this "risk" function, we 
determine the combination of weight percentages for stocks 1 
thru 3 which produces the smallest statistical risk. The above 
described risk can readily be calculated by various means, such 
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as, by a computer program (which is included in Appendix I), 
The output (see Attachment to Appendix I) of the calculation 
includes the weight percent of each stock and the associated 
overall risk level. This calculation can be repeated for a 
range of expected asset return levels and results in generating 
a nonlinear type "bullet" shape defining the limits of minimum 
risk over a range of asset return levels for associated 
standard deviations of funding level (see FIG. 6). The method 
uses historical returns for the plurality of stocks analyzed in 
order to calcuXa^fe the resulting covariance between the 
standard liiability returns and the future returns of the 
potential portfolio of assets. Appendix II illustrates an 
example of a computer program for calculating typical liability 
return data. The method of analysis results in choosing a 
selected set of assets for the portfolio with a strong 
inclination of the selected set of assets to respond in a 
manner such as thie standard asset returns over time, which 
alone can be valuable output. As mentioned hereinbefore, in 
other embodiments, the nonlinear analysis of a plurality of 
assets can involve other methods, such as, index correlation 
parametrization for matching the performance of a target index 
return (see Appendix III) . 

In one embodiment of the invention illustrated in 
FIG. IB, the method is a simplification of the more formal 
procedure of FIG. lA. The method of FIG. IB accomplishes, 
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however, the same result as in FIG. lA but with much greater 
efficiency and speed. 

The method of FIG, IB is broken into eleven steps, and 
the first three steps are substantially the same as the method 
of FIG. lA. The eleven steps and details of each each step are 
described below: 

I- Re^d in Returns of Selected s ecuritie s. Large 
numbers (thousands of securities) can be utilized in this 
procedure. An input file contains the returns for each 
security in the prior periods, such as the previous 24 months. 
Various data is included, such as, security name, industry 
groups code, market capitalization, trading volume, recent 
prices, specific identifiers and estimated bid/ask price 
spreads. The effect of transaction costs associated with the 
spread in the bid/ask price can be included in the performance 
analysis. 

^i- Bead in Returns of Target index tn Trft r k if a 
portfolio is to be constructed for tracking a specific 
financial target index, the returns to that target for the 
relevant period are read into memory arrays. Additional 
identifying information is also read in from the target data 
file. 

Set-UP CPnStrgint S on Upper and I .ower Bounds 
Terms og Per cent of Port f oli o for earh security and each 
Sector pr Jndwptry . a maximum and minimum percentage weight of 
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the portfolio for each security can be specified to constrain 
the portfolio. This can be used to insure portfolio 
diversification an(3[;to control costs associated with trading. 
In addition, sectors of the universe, for example, utility 
stock, can be constrained by maximum and minimum boundaries. 
If one has a single target, one can "short" a stock and take a 
negative minimum position, 

IV. Calculate the Cov ariance of Each Security 

with the Target Index . If a target is used the covariance is 
calculated for each security and stored in an array. If no 
target index is ijji^ed, a zero value for each security is stored 
in the array. 

V. Create an Initial Portfolio bv Selecting 

Highest Covariance Securities and Weighting Them at their Upper 
Limits as Defined bv Selected Constraints. This step creates 
an initial feasible solution to the problem by filling the 
vector of portfolio weights according to the constraints and in 
order of highest covariance. 

VI . . : Calculate the Objec tive Function Value at 
Initial PortfollQ Weight ings. The objective function can be 
defined in a numjtjer of ways. The computer program allows 
monthly or moving quarterly returns to be used for 
optimization. Transaction costs can be considered and their 
importance magnified or reduced relative to other objectives. 
With minor changes, other such goals can be incorporated into 
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the objective function. The key is that once the objective 
function is specified, partial derivatives can be used to guide 
the search for an optimal portfolio. Any example of a 
preference that can be created is an increased weighting for a 
stock with a likely dividend versus one with no dividend, 
VII. Calculate Partial Derivatives for Earh 
Security at Initial Portf olio Weighting s. A partial derivative 
is calculated for each variable (in this case each security), 
and a direction can be determined in which to move the 
individual security weights in order to obtain an improved 
portfolio solution. The partial derivatives are also used to 
determine if the weights are optimal* In the prior art, the 
solution techniques require storage of a full covariance matrix 
array. This storage requirement has limited the practical 
number of securities which could be considered at one time in 
the past methodologies. Indeed, the storage requirements, and 
to a large degree the processing time, varies as the square of 
the number of securities in the portfolio under consideration. 
See the example discussed hereinafter in which the previous 
methodology is compared to the invention. 

In addition, the methods of solution feor these 
problems in the prior art were slow and cumbersome and subject 
to failure when the full covariance array was sensitive or a 
nonunique set of solutions were achievable. This current 
method requires much less storage, uses a rapid solution 
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technique and allows control of the tolerance used for 
optimality. Consequently, the improved methods will select a 
portfolio when several combinations are equally desirable. 

Not only is the amount of necessary memory 
reduced and the computer calculational time greatly reduced, 
there is substantial flexibility in defining the objective, 
assurance of the solution is enhanced and simultaneous 
considerations 6^ large number of securities allows substantial 
improvement in optimising the expected return of the portfolio 
compared to the tiarget< index. 

VIII, Change of Portfolio Weightings bv Moving in 

Direc tion of Improv ement is Indicated bv Partial Derivatives. 
Adjustment of the portfolio weights is achieved by a search 
technique which moves along the constraints and changes in a 
proper direction pfi improvement of the objective function. The 
objective function is calculated at the new weights and a test 
of improvement is made* 

IX, Tegt fog Conveyqence to an OptimpJ, Solution 
Defined by the Kuhn-Tucker Conditions . If the objective 
function is not improving" or if the step size used to adjust 
the portfolio weights becomes extremely small, the search is 
terminated; This solution is normally a Kuhn-Tucker point 
(conventional method of establishing optimality conditions) or 
extremely close thereto within an acceptable epsilon to such a 
point. 
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X. RegYcUnq Cc>n(?i.ti,9ir)g > if the termination 
conditions are not satisfied then one re-calculates the 
objective function value, re-calculates the partial 
derivatives, makes changes in the portfolio weightings to 
> achieve an improved solution and test for convergence. 

The process is by nature iterative and continues 
cycling until a solution is reached in which diminishing 
returns are achieved by further cycling. 

XI. Output of Information. Relevant portfolio 
information is output with security weightings, objective 
function values, purchases and sales necessary to achieve the 
optimum portfolio and industry weightings. 

The detailed output is written into a computer file 
which then can be examined for relevant information. Order to 
buy and sell securities can be developed from the information 
in the output file. 

Included in Appendix IV is an exemplary computer 
software (source output) program illustrating critical steps of 
the method of FIG. IB. Table XI shows exemplary results for a 
program simulation wherein the target index is the Standard and 
Poors 500 stock index. Appendix VI illustrates significant 
distinctions from the optimizer methodology used in the 
copending patent application having serial number 281,560. 



Quantifying Magnitude of Reduced Computation Time 

A test was performed on an IBM compatible PC to 



wo 91/02326 PCr/US90/04328 

^ . . . 

compare the solution speeds of two portfolio optimization 
systems. In a prior system/ the problem of handling large 
numbers of securities in a portfolio selection process 
increased in proportion to the number of securities squared. 
Thus, a problem involving one hundred securities would take 
approximately onQ hundred times as many calculations to solve 
as a problem with ten securities. 

In the current system the solution difficulty 
increases by a factor of less than one times the number of 
securities. In addition, the computing memory required to 
solve the problem is proportional to the number of securities 
rather than, as in the prior system, that number squared. 

Solution Time Goraparisonr 

Hardware: 1MB compatible PC, 386-20Mhz CPU, 80387 
co-prpcessor 

Problem Size: 100 securities 

i 

3% portfolio weight upper bound on each security 
0% portfolio weight lower bound on each security 
Time to solution^ 

Prior system - 9 minutes 34 seconds 

Current system - 52 seconds 

Overall speed-up:. 11.0 x prior system 
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~ 

This speed increase is actually of a larger magnitude 
since approximately 40% of the 52 seconds to solution with the 
current system is spent reading information from external files 
and writing information to other files. The solution time 
speed-up, allowing for reading and writing files, is 
approximately 17.3 x prior system. This speed-up ratio 
increases in proporation to the number of securities considered 
squared. 

The speed increase and the reduced computer memory 
requirements by the current system allow large problems to be 
solved in a short time, requires limited computer memory, and 
uses computer hardware which is relatively inexpensive. 

Circumstances arise regularly in the investment field 
which rapidly change the prospects for securities. The impact 
Of these sudden changes must be incorporated into the security 
valuation system so that rational alternations in the 
investment portfolios may be made. 

Examples of sudden changes include: a company is 
presented with a buy-out offer by another firm; a disaster 
occurs, such as an oil spill, which may impact a firm's stock 
price; monetary or fiscal policy changes are implemented by the 
government, it is important for an investment system to be 
flexible and fast enough to evaluate the impact these changes 
may have on a security portfolio. 
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The current system allows for estimates of partial 
monthly returns to be calculated on any day of the month, for 
these returns to be used in the optimization process, and for 
the results of the analysis to be completed within a few 
minutes. 

One advantage of the current system is that analysis 
of the current investment opportunities can be completed 
rapidly and recommendations for buying and selling securities 
can be quickly generated. This allows investment decisions to 
be made and implemented quickly with confidence. 

Further illustrations of the invention are exemplified 
by various historical simulations shown in FIGS. 2-4 and 
Tables I-III which are taken over the time period of 1975 to 
1987. As listed in Table 1 and in FIGS. 2 and 3, the liability 
stream for a selected pensfion plan can undergo substantial 
variation with time. A portfolio of assets has been analyzed 
in accordance with the preferred statistical method described 
hereinbefore, and detai]^s of the selected portfolio are set 
forth in Tables IV-X. Over the 1975-1987 time period, the 
resulting portfolio of assets shows substantially better 
correlatiop to the liability stream as compared to the Standard 
& Poors 500 return. . Moreover, as seen in Table II and FIG. 4, 
the overall cumuiative return for the portfolio of assets 
selected by the preferred method is far better than the 
Standard & Poors 500. The greatly enhanced stability and good 
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TABLE V 



SWULATION RESULTS 
1975 

HISTORIC SIMUUTION: 
LIABILITY STREAM USED LONG TERM STUDY ABO 



NUMBER 


WEIGHT 


SIC i 


IOC 


SECURITY NAME 


1 


3.00 


67. 


BTC 


BELL CANADA ENTERPRISES 


2 


3.00 


63 


CBB 


CHUBB CORP 


3 


3.00 


33 


X 


USX CORP 


4 


3.00 


35 


HR 


NAVISTAR INTL CORP 


5 


3.00 


54 


WIN 


WINN DIXIE STORES INC 


6 


3.00 


13 


CBI 


CBI INDS INC 


7 


3.00 . 


49 


WUP 


WASHINGTON WTR PWR CO 


8 


3.00 


29 


AHC9 


AHERAOA HESS CORP 


9 


3.00 


35 


UT 


UNITED TELECOMMUNICATIONS 


10 


3.00 


10 


AHX 


AMAX INC 


11 


3.00 


67 


ASA 


ASA LTD 


12 


3.00 


49 


PGN 


PORTLAND GEN CORP 


13 


3.00 


64 


AXO 


ALEXANDER & ALEXANDER SVCS 


14 


3.00 


13 


KMC 


KERR HC6EE CORP 


15 


3.00 


48 


T 


AMERICAN TEL i TELEG CO 


16 


3.00 


63 


CIC 


CONTINENTAL CORP 


17 


3.00 


60 


FBGl 


BANC ONE CORP 


18 


3.00 


33 


NS 


NATIONAL INTERGROUP INC 


19 


3.00 


26 


KHB 


KIHBERLY CLARK CORP 


20 


3.00 


36 


HSU9 


MATSUSHITA ELEC INOL 


21 


3.00 


13 


GAS 


NICOR INC 


22 


3.00 


49 


GPL 


CAROLINA PWR & LT CO 


23 


3.00 


37 


UA 


UNITED TECHNOLOGIES CORP 


24 


3.00 


64 


MMC 


HARSH & MCLENNAN COS INC 


25 


3.00 


12 


BNI 


BURLINGTON NORTHN INC 


26 


3.00 


20 


K 


KELLOGG CO 


27 


3.00 


63 


UFY 


USF&G CORP 


28 


3.00 


34 


AC 


PRIHERICA CORP 


29 


3.00 


13 


SN 


AMOCO CORP 


30 


3.00 


27 


DNY 


DONNELLEY RR & SONS CO 


31 


2.23 


10 


AL 


ALCAN ALUM LTD 


32 


• 1.92 


33 


IAD 


INLAND STL INDS INC 


33 


1.63 


32 


PPG 


PPG INDS INC 


34 


1.41 


60 


FHL 


FIRST BK SYS INC 


35 


1. 00 


49 


PLT 


PACIFIC LTG CORP 


36 


0.83 


67 


LEH 


LEHMAN CORP 


37 


0.51 


10 


UNP 


UNION PAC CORP 


38 


0.48 


13 


OXY 


OCCIDENTAL PETE CGPCR 
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SIMULATION RESULTS 
1980 

HISTORIC SIHULATIOH: 



1 

2 
3 
4 
5 
6 
7 
8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 



LIABILITY STREAM USED 



NUMBER WEIGHT 



3.00 
3.00 



00 
00 
00 
00 

oio 

PD 
3.00 
3.00 
3.00 
3.00 
3. CO 
3.00 
3.00 
3.0b 

3. bo 

3.00 
3.00 

3.Qa. 

3.00 

3.00 

3.00 

3,b0 

3.00 

3.00 

3:00 

3;oo 
3;bo 

2.33 
2.06 
2.05 
liOO 
1.00 
0.95 
0.62 



cfr £ 


IDC 




nriN 


28 . 


Irr 


60 ' ' 




36 


UUP 


60 


nUD 


25 




28 


WLn 


20 


DAI 
KnL 


49 




28 


DAD 


A^ 

26 


lirr 


r A 

60 




A A 

20 




30 


rT 


49 


DOM 

rUn 


V A 

13 


1 1 Y 


Zo 








33 




26 


rnr 


A A 

20 


DM 


20 






AC 


28 


GLXO 


34 


GS 


37 


GH 


23 


ISS 


48 


T 


20 


CPB 


36 


GNB 


37 


HL 


37 


F 


36 


C 


49 


PPL 


60 


AXP 


32 


OCF 


20 


PSY 



FR 



LONG TERM ABO 

SECURITY NAME 

NIAGARA HOHAVnC PMR CORP 
INTERNATIONAL FLAVORS J 
HELLOM BANK CORPORATION 
WHIRLPOOL CORP 
NORHEST CORP 
KIHBERLY CLARK CORP 
WARNER LAMBERT CO 
RALSTON PURINA CO 
KANSAS GAS i ELEC CO 
BORG WARNER CORP 
UNION CAMP CORP 
CHEMICAL HEMYORK CORP 
KELLOGG CO 

GOODYEAR TIRE & RUBR CO 
POTOMAC ELEC PMR CO 
LOUISIANA LO i EXPL CO 
PROCTER i GAMBLE CO 
FEDERATED DEPT STORES INC 
NATIONAL INTERGROUP INC 
FORT HOWARD CORP 
BORDEN INC 
GENERAL HLS INC 
PRIHERICA CORP 
GLAXO HLDGS PLC 
GILLEHE CO 
GENERAL HTRS CORP 
INTCRCO INC 

AMERICAN TEL i TELEG CO 
CAMPBELL SOUP CO 
GOULD INC 

MARTIN HARIEHA CORP 
FORD HTR CO DEL 
CHRYSLER HLDG CO 
PENNSYLVANIA PWR & LT CO 
AMERICAN EXPRESS CO 
OWENS CORNING FIBERGLAS CO 
PILLSBURY CO 
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Code 

10 
12 
13 
20 
26 
27 
29 
32 
33 
34 
35 
36 
37 
48 
49 
54 
60 
63 
64 
67 



1975 INDUSTRY WEIGHTINGS 

Electric, Gas and Sanitation Utilities, and 
Banking Restricted to 10% of Portfolio. 



Industry 
Metal Mining 

Bituminous Coal & Lignite Hin. 
Oil i Gas Extraction 
Food Kindred Products 
Paper and Allied Products 
Printing Publishing and A.P. 
Petroleum Refin. i Rel. Prod. 
Stone Clay Glass i Cone. Prod. 
Primary Metal Industries 
Fab. Metal Prod. Ex. M.iT.E. 
Machinery Except Electrical 
Ele. and Ele. Hach. 
Transportations Equipment 
Communication 

Electric Gas and Sanlt. Serv. 

Food Stores 

Banldng 

Insurance Carriers 
Ins. Agents Brolc. Serv. 
Hold, and Other Inv. Comp. 



% of 
Portfolio 

5.74% 
3.00 

12.48 
3.00 
3.00 
3.00 
3.00 
1.63 
7.92 
3.00 
6.00 
3.00 
3.00 
3.00 

10.00 
3.00 
4.41 
9.00 
6.00 
6.83 



SPANNING PORTFOLIO 
1/1/75 



Market Capitalization 
Average 
High 
Low 

Shares Outstanding 
Average 

Dividend Yield 

Price Earning Ratio 

Number of Stocl(S 

Turnover 1975-1976 



(000,000) 
SI. 203 
$8,393 
S 137 

(000) 
74,840 shares 

5. IX 

13. Ix 

38 

21. 3X 
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Code 

13 
20 
23 
26 
28 
30 
32 
33 
34 
36 
37 
48 
49 
53 
60 



1980 INDUSTRY WEIGHTINGS 

Electric. Gas and Sanitation Utilities, and 
Banking Restricted to 10% of Portfolio. 



Industry 

Oil , , 

Food Kindred Products 
App. i 0th. Fin. Prod. 
Paper and Allied Products 
Chemical and Allied Prod. 
Rubber and Hisc. Plast. Prod. 
Stone Clay Glass i Cone. Prod. 
Primary Metal Industries 
Fab. Metal Prod. Ex. M.iT.t. 
Ele^ and Ele, Mach. 
Transportations Equipment 
Conmunication 

Electric Gas and Sanit. Serv. 
General Herch. Stores 
Banking 



% of 
Portfolio 

3.00% 

15.62 
3.00 
9.00 

15.00 
3.00 
0.95 
3.00 
6.00 
8.05 
7.38 
3.00 

10.00 
3.00 

10.00 



SPANNING PORTFOLIO 
1/1/80 



Market Capitalizatton 
Average 
High 
Low 

Shares Outstanding 
Average 

Dividend Yield ]. 

Price Earning Ratio 

Number of Stocks 

Turnover 1980-1981 



(000,000) 
$1,790 

513,311 
S 247 

(000) 

130,229 shares 
6.9% 
6.6x 
37 
30.4% 
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TABLE IX 

SIMULATION RESULTS . 
1985 

HISTORIC SIMULATION: 



LIABILITY STREAM USED LONG TERM STUDY ABO 



NUMBER 


WEIGHT 


SIC i 


luc 




1 


3.00 


bo 


ITD 


PPMNFY ^ r INC 


2 


^ Art 
3 .00 




UP 


HAUATTAN ELEC I 


3 


3.00 


CD 

bo 


urn 




4 


3.00 






^nilTHFAST BKG C 


5 


^ AA 

3.00 




DMf- 


ixUOdCLL i#wrvr 


0 


^ AA 

3*00 


Zo 


A\/D 
AYr 


AvriN ppon^ INC 


7 


^ A A 

3.00 


60 


nr 1 


MANIIPArTlIPrP^ RAN 


8 


3.00 


AA 


TdILl 


LAon. IDlLLo U.V 


9 


A A 

3.00 


A A 

49 


not 

DPL 


noi Twr 
UrU INt 


10 


A A 

3.00 


3b 


net 


MTi 1 Tpnor ropp 


11 


3*00 


A A 

20 


Iff 
UL 


fiifTft/ro Pt r 


12 


<9 A A 

3.00 


3o 






13 


3.00 


Z7 




rnMMPorr n fapt 


14 


3.00 


28 


CAD 


cniiTRP mpD 


15 


3.00 


1 A 

28 


PFE 


DCTTCD TMH 

PrlitK INU 


16 


3.00 


A A 

20 


t^A 

KO 


LOLA tULA tU 


17 


3.00 


45 


A A 1 

DAL 


nri TA ATO 1 Turc 
UcLTA AiK Linbb 


18 


3.00 


67 


NES 


ftfcii ruAt Atin n c 
NEW ENGLAND lLl 


19 


3.00 


AA 

89 


CI t 
SU 


STONc a VftDbttK 


20 


3.00 


61 


Al tC 

GwF 


GREAT WtSTN rIN 


21 


3.00 


28 


SKL 


SnlTHKLlNE DttNrlAn 


22 


3.00 


28 


SRG 


SCHERING PLOUuH 


23 


3.00 


51 


SYY 


SYSCO CORP 


24 


3.00 


38 


IIUU 

nnn 


MINNESOTA HNu & 


25 


2.93 


78 


A V A 

DIS 


nTPfcfPV UAI T AA 

DISNEY HALT CO 


26 


2.68 


25 


■ tA 

HB 


HILLENBRAND IND 


27 


2.66 


49 


B6E 


BALTIHORE GAS & ELE 


28 


2.53 


38 


KDO 


HEDTRONIC INC 


29 


2.30 


73 


t/ A*. 

K6S 


KELLY SVCS INC 


40 






PPMP 


PPTMP COMPUTER 


31 


1.77 


28 


SYN 


SYNTEX CORP 


32 


1.76 


16 


. FLR 


FLUOR CORP 


33 


1.70 


60 


FBGl 


BANC ONE CORP 


34 


1.54 


63 


AHN 


AHHANSON H F & 


35 


1.35 


49 


AYP 


ALLEGHENY PUR SYS 


36 


1.22 


48 


CTX 


CENTEL CORP 


37 


1.09 


20 


CFG 


CPC INTL INC 


38 


0.77 


28 


PG 


PROCTER & GAMBLE 


39 


0.59 


28 


AHP 


AMERICAN HOHE PRODS 


40 


0.54 


34 


GS 


GILLEHE CO 


41 


0.41 


53 


DNC 


DAYTON HUDSON C 


42 


0.13 


47 


FEDX 


FEDERAL EXPRESS 
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Code 

16 

20 

23 

25 

27 

28 

34 

35 

38 

45 

47 

48 

49 

51 

S3 

58 

60 

61 

63 

67 

73 

78 

89 

99 



3/ 

TABLE X 

^ 1985 INDUSTRY WEIGHTINGS 

Electric, Gas and Sanitation Utilities, and 
eSinJ Restricted to lOX of Portfolio. 

Industry 

Const.. 0th Than B.C.-G.C. 
Fbod-Wndred Products 
i Aip. i 0th. Fin. Prod. 
Farniture and Fixtures 
teinting Publishing and A. P. 
Oiemical and Allied Prod. 
Fab. Hetal Prod. Ex. H.iT.E. 
MaQhinery Except Electrical 
Heasv Anal, i Cent. Inst. Etc. 
i.: Transportation By Air 
■<{• Tranipbrtation Services 
""• Cotaptini cation 

Electric Gas and Sanit. Serv. 
tftolesale Trade-llondur. Goods 
General Herch. Stores 
Eating and Drinking Places 

Iredl^en. 0th. Than Banks 

Insurance Carriers 

Hold. And Other Inv. Coop. 

Business Services 

Motion Pictures , 

Miscellaneous Services 

HonclassifiaMe Establishments 
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% of 
Portfolio 

1.76% 
7.09 
3.00 
2.68 
3.00 
18.13 
0.54 
5.05 
8.53 
3.00 
0.13 
1.22 
10.00 
3.00 
3.41 
3.00 
7.70 
3.00 
1.54 
3.00 
2.30 
2.93 
3.00 
3.00 



SPANNING PORTFOLIO 
1/1/85 



Market Capitalization 
Averigft 
High 
Low 

Shares Outstanding 
Average 

Dividend Yield 

Price Earning Ratio 

Number of Stocks 

Turnover 1985-1986 



(000.000) 
$3,045 
S11.689 
J 373 

(000) 
97.120 shares 

4.0% 

ll.Ox 

41 

15.0% 
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Table XI 
ANALYSIS OF PROGRAM SIMULATION 

Target & Index represent the Standard & Poors 500 Stock Index 
portfolio represent optimizer chosen portfolio. 
One example of the benefits of the technique are seen in the 
monthly statistics, where the annualized standard deviation of 
returns is 14.6% for the optimized portfolio vs. 16.04 for the 
S&P 500. 



ANNUALIZED MONTHLY STATISTICS 



Target 
Portfolio 
Portfolio Diff 
Index 

Index Diff 



MEAN 
15.88% 
17.97% 
2.09% 
15.88% 
.00% 



STD DEV 
16.04% 
14.60% 
4.40% 
16.04% 
.00% 



SKEW * lo'^e 
-50.00 
-19.72 

.01 
-50.00 

.00 



ANNUALIZED QUARTERLY STATISTICS 



Target 
Portfolio 
Portfolio Diff 
Index 

Index Diff 
ANNUAL STATISTICS 



MEAN 
16.18% 
18.40% 
2.23% 
16.18% 
.00% 



STD DEV 
17.04% 
16.19% 
3.99% 
17.04% 
.00% 



SKEW * lo'^e 
-172.12 
-41.29 
1.58 
-172.12 
.00 



Target 
Portfolio 
Portfolio Diff 
Index 

Index Diff 



MEAN 
16.38% 
19.01% 
2.63% 
16.38% 
.00% 



STD DEV 
13.63% 
13.43% 
4.43% 
13.63% 
.00% 



SKEW * lo'^e 
-670.98 
-61.67 
-23.31 
-670.98 
.00 
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statistical correlation with the liability return is further 
evident in Table III ^and FIG. 5, wherein detailed comparisons 
are made between the selected portfolio of assets and the 
standard liability return. 

In fig!" 6 a range of simulation funding returns for 
the portfolio of assets are compared with a typical pension 
fund a mixture of stocks, bonds, real estate and treasury 
bills. Clearly, the risk is much higher for the typical 
pension fund; and dramatic improvement in the return, or 
reduction of risk, results when only 35% of the standard 
pension fund is modified using the method of the invention. 

In another embodiment, a portfolio of assets can be 
constructed by selecting a portion of a total portfolio with 
assets having optimal correlation of asset return to a 
liability or finaincial index. The remainder of the portfolio 
comprises futures contracts which are combined with the 
correlated portfolio portion to achieve in effect an optimum 
correlation for the entire portfolio of assets. Further 
details are set forth in Appendix IV. 

In another aspect of one embodiment, control can be 
exerted over pension plan surplus by adjusting the level of 
risk selected for a portfolio of assets. As illustrated in 
FIG. 6, the expected return can be selected at various levels 
with the degree of risk, or standard deviation of the funding 
level, generally increasing as one moves from a position of 
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minimum risk at the top of the "bullet" to higher future 
returns. Control over a pension plan surplus, or for that 
matter any plan for which you wish to respond dynamically to 
control risk/return in concert, can be accomplished over a wide 
range of risk and return values. Such an approach can be used 
to manage return under variable risk positions and minimize 
insurance costs for protecting against underfunding of a plan, 
such as falling below a predetermined minimum floor. 
Consequently, as the funding level approaches 100% a minimum 
risk portfolio of assets should be constructed using the 
methods described hereinbefore. As the surplus accumulates, 
the acceptable risk level can be increased for the portfolio of 
assets by dynamic modification of the portfolio asset 
components. One can utilize futures contracts as an overlay 
for an underlying portfolio of assets, having been selected by 
the basic invention described previously, to create in effect 
an optimum statistical correlation for the entire portfolio, 
including the futures contracts. As the surplus approaches 
10% - 20% excess, a portfolio of assets can be constructed 
resulting in a much higher level of future return. For 
example, in FIG. 6, the change in future return from minimum 
risk to the highest return data point is about a 35% greater 
return but with an accompanying 70% - 80% increase in standard 
deviation compared to the minimum risk point. 
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In a further embodiment of the invention one can 
overcome problems associated with an indefinite covariant 
matrix. Th^ current system and method does not use the full 
covariance matrix bO: extract a partial derivative to guide the 
search process, and thus it is not subject to one of the 
failures that results using standard quadratic programming. 

To us standard quadratic programming algorithums the 
covariance raatrizi'must be positive semi-definite, or positive 
definite. This means technically that no row of the matrix can 
be replicated by a linear combination of other rows. 

However, this condition occurs when there are fewer 
returns than the number of securities under consideration. 
This is a significant shortcoming of the standard methodlogy. 
For example, to consider five hundred securities 
simultaneously, the user must supply at least five hundred and 
one returns for each security being considered. If monthly 
data is being use^ for the return series, at least forty-one 
years of data musC/'be available for each security under 
consideration r Most securities have not been in existence for 
this period of tiipe. 

If not enough returns are available or if the matrix 
is indefinite, the standard quadratic solution techniques will 
fail to find a solution point that is optimal. 

The current system allows a solution, which is at 
least as good as any other feasible solution, to be achieved 
without regard to the sensitivity of the covariance matrix. 
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Problems that otherwise could not be solved by 
standard technologies are solved by the current system. This 
allows practical portfolios to be selected even though there is 
limited available return information. 

In another aspect of the invention involving practical 
applications of the current methodology to managing security 
portfolios, the ability to invest dividend income and other 
cash flows efficiently is also an important element in 
effective management. The current system allows a portfolio 
manager to invest available cash in the most efficient 
securities while considering the current portfolio holdings. 
Effectively this allows the manager to invest in securities 
which best enhance the current portfolio position without 
selling any of the current holdings. In a practical portfolio 
management system the ability to reinvest cash flows 
efficiently is always an important consideration. The current 
system provides this ability. An illustration of a 
reinvestment solution is provided in Appendix V which lists 
data used and results obtained in performing the reinvestment 
method. 

In addition to finding optimal portfolios for tracking 
financial targets, the current system allows other objectives 
to be considered and incorporated in the solution. Examples 
include supplementing the basic objective function with an 
income objective, tax impact objective, or a company cash flow 
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objective. The system is flexible enough to allow the 
objective function to be customized for particular applications. 

As an example, an investor who has a preference for 
securities with high dividend yields will specify an objective 
function which explicitly states the trade off between 
portfolio tracking and dividend income. The computer routine 
to optimize a portfolio (the optimizer) will extract the 
partial derivatives of this objective function and proceed to 
select an optimal portfolio which exactly incorporates the 
specified trade-off between dividend income and tracking. This 
investor then has a custom solution which addresses his 
particular concerns and requirements. 

The consideration of transaction costs associated with 
buying and selling securities is incorporated into the 
portfolio optimization system to control expenses due to 
trading. The trade off between tracking accuracy and 
transaction costs can be specified by the investor. This 
allows for a customized objective function, with regard to 
expected transaction costS/ for each client. 

Investors may also have tax effects to consider when 
trading from one security to another. The taxable gains and 
losses and their impact on expected return can be specified by 
the investor so. as to control these costs. 

Another preference which can be incorporated into the 
objective function and handled explicitly by the optimization 
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system is the consideration of cash flows. An investor may 
have a preference for investing in companies which have large 
and positive cash flows. This objective can be incorporated 
into the system and resulting portfolios will reflect this 
investment goal. 

While preferred embodiments of the present invention 
have been illustrated and described, it will be understood that 
changes and modifications can be made therein without departing 
from the invention to its broader aspects. Various features of 
the invention are defined in the following claims. 
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APPENDICES TO SPECIFICATION 



APPENDIX I - 



Spanning Program for creating an optimal 
target tracking portfolio of securities 
(Appendix I ppV 1-271). Also includes 4 
pages of Output. 



APPENDIX 11^- 



Li ability Return Program for creating a 
return series for a liability stream from 
yields (Appendix II pp, 1-g). 



APPENDIX III - Correlation Portfolio Program for creating 
an optimal index correlation portfolio with 
securities (Appendix III pp. m^o)'. 



APPENDIX IV - 



APPENDIX V - 



"Fasttrack" Program for analyzing large 
numbers of securities in a rapid, efficient 
manner to provide optimum correlation of 
asset return to a time dependent financial 
index. Substantial computer memory storage 
reduction is also achieved. 

Data listing of method for reinvestment of 
available cash flow beginning from a current 
portfolio. 



APPENDIX VI - 



Example illustrating distinctions between 
Current system and prior art. 
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APPENDIX I 

SPANNING PROGRAM 
To create an optimal target tracking portfolio of securities. 

(Appendix I pp, 1- 27) 



1 
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C UPDATE December 6, 1988 

Q 

PROGRAM SPAN 

C **^ 

INCLUDE 'COMMON. F' 
CHARACTER*30 HOLD 
LOGICAL I COUNT 
LOGICAL HOLOEX 



C ************* WRITE COPYRIGHT TO SCREEN 
INCLUDE 'COPYRIGHT. F' 

C 



**************** 



******** 



C ************* Qi^jp^ j^py.^ SECTION ************************ 

5 CONTINUE 



C 



DEFAULT READ IN OF 24 MONTHS PRIOR = NSTAT. NRESTS FOR B( ) 



**** 



NSTAT " 24 
NRETS = 24 
BULLET =0.0 
TARGET = 1. 
TURN = 100. 
ITYPE = 1 
ICHANGE = 0 
XBLOWER = .ODO 
XBUPPER = 3.0D0 
SMAX = XBUPPER 



C *7rrnAT******* CALL FLASH TO INTRODUCE SPAN **** 



CALL FLASH 
************* 



**** 



r **1L1P"^'^ -^QV. .FALSE.) THEN 

r«i . T.rrn« ^'^^^ ^^TRO TO INPUT PARAMETERS 

C CALL MODIFY TO CHANGE PARAMETERS 

' *"lFaH*;~^ READ IN THE HOLDINGS FILE 

" READ(i>nPOrHOLD"'''°'-^° *' ^"^^^^^^^ "^E' 

INQUIRE (FILE = HOLD. EXIST = HOLDEX) 



**** 
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vwvwwirwii 



c *** 



* 
* 



C *** ND 

C *** ID 

C *** N2D 

C *** KEQ 



C *** 
c *** 



NRT 
NST 



TOTAL NUMBER OF STOCK WEIGHTS BEING COMPUTED 
THE MAXIMUM NUMBER OF CONSTRAINT EQUATIONS 
THE DIMENSION ON WORKING AREAS ( 2 * NO ) 
NUMBER OF CONSTRAINTS WHICH WILL BE EQUALITIES 

NUMBER OF RETURNS 
USED FOR STATISTICS 



C NDID^^ iJoVli ' ' ^Hl^^Zm 46 

C ********* pQ,^ ^f^i^^ Jjy^j .^^ 

c *** 

C LAST UPDATED: November 19, 1988 
IMPLICIT REAL (A-H,0-Z) 
IMPLICIT INTEGER (I-N) 
PARAMETER (KEQ - 0) 
280 ^ ^ I 

PARAMETER (ND - 280, ID - 6 
PARAMETER (NDID - 286, ISIZEQ 

PARAMETER (ND « 250, ID - 6 
PARAMETER (NDID - 256, ISIZEQ 

100 

PARAMETER (ND - 100, ID - 6 
PARAMETER (NDID - 106, ISIZEQ 



* 

* 
* 



N2D - 560} 
- 645246) 

N2D - 500) 
' 516126) 



N2D - 200) 
86526) 



PARAMETER (NRT - 48, NST - 3) 
PARAMETER (LENGTH - 40) 

character*(length) fileout,title,name{nd),portrets garb answer 
& old.lianame.liabin.stockinIKdinjocw^^^^ 

common /arry/ a(nd,nd),c(nd,io) 

COMMnil ^bctK ite^iSI •!P'-<NDi.BDU(ND),X(ND).OLDWT(ND) 

COMMON /TRA2/ PRICE(ND),VOL(ND) CUM(ND) VARff)^^ 
COMMON /RTNIN/FILEOUT.L ABIN.St6H;KN l6c 
^COMMON /DAT/ NMSAVE.^W.NSIMlNSTOCKsSra^ 

COMMON /MAX/ ST0CKMIN,YIMAX,SHAX.SCALE.BRET XFACTOR XBIIPPFR 
C END OF co«SiT"''"''"'''*="'™'^''^"-'«"'^^^ 
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APPENDIX II 

LIABILITY RETURN PROGRAM 
To create a return series for a liability stream from yields. 

(Appendix II pp. 1-2) 
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APPENDIX III 

CORRELATION PORTFOLIO PROGRAM 
To create an optimal index correlation portfolio with securities. 

(Appendix III pp. 1-30) 
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C *** PARAMETERS ARE USED TO DIMENSION THE COMMON ARRAYS * 
C **♦ IN EACH SUBROUTINE. * 
C *** * 

C *** ND « TOTAL NUMBER OF STOCK WEIGHTS BEING COMPUTED. * 

C *** ID - THE MAXIMUM NUMBER OF CONSTRAINT EQUATIONS. * 

C *** N2D - THE DIMENSION ON WORKING AREAS ( 2 * ND ) * 

C *** KEQ - NUMBER OF CONSTRAINTS WHICH WILL BE EQUALITIES * 

C *** NRT - NUMBER OF RETURNS * 
C *** NST - USED FOR STATISTICS * 

C *** ISIZEQ • 2 * ( ND*2 + ID )**2 + 8*{ND*2+ID) +6 * 
C *** NDID - ND + ID * 

C LAST UPDATED: November 19, 1988 
IMPLICIT REAL (A-H,0-Z) 
IMPLICIT INTEGER (I-N) 
PARAMETER (KEQ - 0) 

C 260 

C PARAMETER (ND " 280, ID > 6 , N2D - 560) 
C PARAMETER (NDID - 286, ISIZEQ - 645246) 
C 250 

C PARAMETER (ND - 250, ID - 6 , N2D - 500) 
C PARAMETER (NDID - 256, ISIZEQ ' 516126) 
C 100 

' PARAMETER (ND - 100, ID - 6 , N2D - 200) 
PARAMETER (NDID • 106, ISIZEQ - 86526) 

PARAMETER (NRT - 48, NST - 3) 



PARAMETER 



LENGTH - 40) 



'1 

CHARACTER*(LENGTH) FILEOUT,TITLE,NAME(ND),PORTRETS,GARB,ANSWER, 
& 0LD,LIANANE,LIABIN,ST0CKIN,B0NDIN,IDC(ND),SIC(99) 



COMMON /km/ A(ND,ND),C(ND,ID), H(ND,NO), G(ND), F(ND) 
COMMON /VECT/ B(ND),D(ID),BDL(ND),BDU(ND ,X(ND),0LDWT(ND) 
COMMON /RETS/ RET(NRT,ND),ANIH(NRT),PRET(NRT),!SIC(NO),AVE(ND) 
COMMON /TRAC/ COV NST,NST),STAT(NRT,NST),COREL(NRT,NST) 
COMMON /TRA2/ PRICE(ND),V0L(ND),CUM(ND),VAR(3),LIANAHE,TITLE 
COMMON /RTNIN/FILEOUT,LlABIN,STOCKIN,BONDIN,IDC,NAHE,SIC 
COMMON /DAT/ NMSAVE,NM,NSIHS,NSTOCKS,NBONDS,NSTAT,NRETS,IDUHMY, 
& IHOLD 

COMMON /MAX/ STOCKHIN,YIMAX,SMAX,SCALE,BRET,XFACTOR,XBUPPER, 
& XBL0WER,BULLET,TAR6ET,TURN, BULLl ,TARG1,TURN1 .PORTVAL 
: END OF COMMON. F 
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APPENDIX IV 



-FASTTRACK- PROGRAM FOR ANALYZING 
LARGE NUMBERS OF SECURITIES IN A RAPID, 
EFFICIENT MANNER TO PROVIDE OPTIMUM 
CORRELATION OF ASSET RETURN TO A TIME DEPENDENT 
FINANCIAL INDEX. SUBSTANTIAL COMPUTER 
MEMORY STORAGE REDUCTION IS ALSO ACHIEVED. 
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ENDIF 
CONTINUE 
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WRITE(*,*)' ' WARNING OLDWTS TOTAL < .99, CONTINUING ' 

WRITE(*,*)' OLDTOT - '.OLOTOT 
ELSEIF(0LDT0T.GT.I.01)THEN .„...r»...T„^ / 

WRITE(*,*)' WARNING OLDWTS TOTAL > 1.01, CONTINUING 

WRITE(*,*)' OLDTOT - '.OLOTOT 
ENDIF 

NGONE-NCOUNT 

IF (IDCTEST.EQ.l) THEN ^ 

WRITEf*,*)' A TOTAL OF ',VALTURN*100.0, ' % OF PORT SOLD' 
WRITE{*,*)NGONE,'SYr«.S NOT FOUND, PROGRAM CONTINUING' 

ENDIF 

r . IF THE 1ST TIME THROUGHT SET OLDWT()-0 

ELSEIF(ISKIP.EQ.1}TKEN 

DO 98 I-1,NUM 
OLDWT{I)-0.0 

CONTINUE 
ENDIF 

******** IF MULTI RUN (iOLO) READ WTS FRO OLD. 

ELSEIF(IFLAG{l):Ei|.3)THEN " i 

0PEN(21,FILE - 'OLD', STATUS - 'OLD') 
DO 94 I - l.NUH 

READ{Zl,*)OLDWT(i) 

IF (TURN.EQ. Cr.0)OLOWT(I)-0.0 

THIS NEXT LINE ADDED 4/18/69 
OLDWT{I)-0.0 

OLDWT(I) - OLDWT(I)/100. 
CONTINUE 
CL0SE(21) 
ENDIF 

RETURN 

END 

END or SUBROUTINE HOLDING ************** 
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APPENDIX V 



REINVESTMENT OF AVAILABLE CASH 
BEGINNING FROM CURRENT 
PORTFOLIO HOLDINGS 
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29 PETROLEUM REFIN. & REL. PROD. 3.21 % 

30 RUBBER AND MISC. PLAST. PROD. 2.70% 
32 STONE CUY GLASS & CQNC. PROD. 0.45 % 

35 MACHINERY EXCEPTILECTRICAL 2.09 % 

36 ELE. AND ELE. HACH. 3.00% 

37 TRANSPORTATIONS EQUIPMENT . 2.20 % 

38 HEAS. ANAL. & CQSIT. INST. ETC. 6.45 % 
45 TRANSPORTATION BY AIR 3.21 % 

48 COMPIUNICATION 1.69 J 

49 ELECTRIC GAS AND SANH. SERV. 17,31 % 
52 BUILDING HAT. HARD. GAR.' SUPP. 3.17% 

60 BANKING , t 1-56 * 

61 CREO. AGEN. OTH. THAN BANKS. 2.83% 
65 REAL ESTATE - 2.69 % 
67 HOLD. AND OTHER INV. COHP. 17.49 % 
73 BUSINESS SERVICES 4.07% 
75 AUTO. REPAIR SERV. AND GAR. 1.19 % 
80 HEALTH SERVICES 3.63% 

TOTAL - 99.99 % 
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This result Is identically the same as extracting the partial derivative 
from the full covariance matrix. However, fewer calculations are required 
and less computer memory is needed to obtain the partial derivatives. 
Mathematically, the results proy^ to be exactly the same. 

Once the partial derivatives have been extracted for each of the 
variables, it is a siftiple matter to search for a better solution using the 
partials to indicate a direction which Improves the objective function. 
The simple process of letting the partial derivatives guide a search, 
recalculating the partials at the new solution, and searching again, leads 
quickly to an optimal solution. The search proceeds downhill until the 
objective function cir\ no longer improve. This is the optimal solution. 

In summary, the current system requires fewer calculations to determine 
an optimal point and therefore can solves problems in much less time. The 
current system alsa.requi res much less computer memory and therefore can solve 
much larger problem. % analyzing more securities at once, a better solution 
can be achieved. 
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OUTPUT FROM PROGRAM 
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1. A method of correlating a future asset return of a 
portfolio to future financial liabilities variable over time, 
comprising the steps of: 

(a) determinihg future payments needed over time to 
fulfill said future financial liabilities; 

(b) selecting a standard index having a standard 
asset return oviBt time for meeting said future payments; and 

(c) analyzing a plurality of assets for identifying 
selected ones of said assets for providing optimum correlation 
of the future return of said selected assets to said standard 
asset return; said optimum correlation achieved by calculating 
a minimum standard deviation for the difference between the 
future returns of said plurality of selected assets and said 
standard asset return, ranking said selected assets from 
smallest to largest covariance, making an iterative change in 
weights of said selected assets, reordering the rank of said 
selected assets and continuing said iterative weight and rank 
change of said selected assets until the covariance converges 
to an optimal solution^ 

2. The method as defined in Claim 1 wherein said future 
financial liabilities comprises a pension plan. 



wo 91/02326 PCr/US90/04328 



3. The method as defined in Claim 1 wherein said plurality of 
assets exhibit sensitivity to the same parameter as said future 
financial liabilities for which said time dependent optimum 
correlation is determined. 

4. A system for correlating a future asset return of a 
portfolio to future financial liabilities variable over time, 
comprising : 

means for determining future payments needed over time 

to fulfill said future financial liabilities; 

means for selecting a standard having a standard asset 

return over time to meet said future payments; and 

means for analyzing a plurality of assets for 

identifying selected ones of said assets for providing optimum 

correlation of the future return of said selected assets to 
said standard asset return, said optimum correlation achieved 
by means for calculating a minimum standard deviation for the 
difference between the future returns of said plurality of 
assets and said selected standard asset return, said 
calculating means further ranking said selected assets from 
smallest to largest covariance, making an iterative change in 
weights of said selected assets, reordering the rank of said 
selected assets and continuing said iterative weight and rank 
change of said selected assets until the covariance converges. 
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5. The system as defined in Claim 4 wherein said future 
financial liabilities comprise insurance liabilities and 
pension plan liabilities. 

6. The method as defined in Claim 4 wherein a computer 
performs the functional operations. 
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1/7 - 

PREPARE INPUT: CREATE FILES WHICH CONTAIN RETURN INFORMATION 
ON THE SECURITIES TO BE USED. THE TARGET RETURNS. SECURITY 
IDENTIFIERS. CURRENT PRICE. MARKET CAPITALIZATION AND CURRENT 
PORTFOLIO WEIGHTS.' "ESTABLISH STATISTICAL PROPERTIES 
OF LIABILITIES." 



I 



SET UP PROBLEM: SPECIFY NUMBER AND TYPE OF SECURITIES, TIME 
PERIOD, NUMBER OF RETURNS USED TO CALCULATE COVARIANCE AND TO 
CALCULATE AVERAGES, DISTANCE TO MOVE UP AND DOWN COVARIANCE BULLET. 
SCALING FACTOR FOR TARGET. AND FACTOR FOR CONTROLLING TURNOVER. 



I 



.Jl^L^l^^^^^^^ INFORMATION: SECURITY NAMES. IDENTIFIERS, 
INDUSTRY CODES, PRICES, MARKET CAPITALIZATION, WEIGHT IN CURRENT 
PORTFOLIO, PERIOD RETURNS; TARGET NAMES, TARGET PERIOD RETURNS. 



I 



CALCULATE SPANNING COVARIANCE ARRAY: DEFINED AS THE 
COVARIANCE OF SECURITY RETURN LESS TARGET RETURN WITH ALL 
OTHER SECURITIES RETURNS LESS TARGET RETURNS 
IN SPECIFIED RETURN PERIOD. 



I 



CALCULATE AVERAGE RETURN FOR EACH SECURITY. 



I 



SET UP CONSTRAINTS ON SECURITY TYPE, SECTOR 
TYPE, INDIVIDUAL SECURITY WEIGHTS. 



FORMAT PROBLEM IN STANDARD FORM FOR SOLUTION BY 
GENERALIZED QUADRATIC PROGRAMMING TECHNIQUE. 



I 



^^^"^L ^P^^"^" PROBLEM WITH QUADRATIC PROGRAMMING 
SOLUTION ALGORITHM. CHECK SOLUTION FOR FEASIBILITY. 



I 



crfiK^o^t'^IL^^^^^S^^^S: COVARIANCE OF PORTFOLIO, ALPHA. BETA, AND 
STANDARD ERROR WITH TARGET IN SIMULATION PERIOD; PORTFOLIO RETURNS 
°^SSL''^''^°''• PORTFOLIO RETURNS IN SIMULATION PERIOD, 
TURNOVER, SECTOR WEIGHTS, AND SECURITY TYPE WEIGHTS. 



I 



rSSSlH? SECURITY WEIGHTS, SENSITIVITIES. 
IDENTIFIERS. AND NAUFS: <;TATT<;TTr<: tmpiit PAOiiirTroc 

SORTED BUYS AND SELLS, ANO" SECTOR WEIGHTS" 



SUBSTITUTP cucrrr 
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PREPARE INPUT: CREATE FILES WHICH CONTAIN RETURN INFORMATION ON 
THE SECURITIES TO BE USED. THE TARGET RETURNS SECURITY IDENTIFIERS, 
CURRENT PRICE, MARKET" CAPITALIZATION AND CURRENT PORTFOLIO 
WEIGHTS. "ESTABLISH STATISTICAL PROPERTIES OF LIABILITIES." 



I 



SET UP PROBLEM: SPECIFY NUMBER AND TYPE OF SECURITIES, TIME 
PERIOD, NUMBER OF RETURNS USED TO CALCULATE COVARIANCE AND TO 
CALCULATE AVERAGES, DISTANCE TO MOVE UP AND DOWN COVARIANCE BULLET, 
SCALING FACTOR FOR TARGET, AND FACTOR FOR CONTROLLING TURNOVER. 



I 



READ IN NECESSARY INFORMATION: SECURITY NAMES, IDENTIFIERS, 
INDUSTRY CODES, PRICES. MARKET CAPITALIZATION. WEIGHT IN CURRENT 
PORTFOLIO. PERIOD RETURNS: TARGET NAME. TARGET PERIOD RETURNS. 



I 



DEFINE VARIANCE AS SUM OF SQUARED DIFFERENCES BETWEEN 
PORTFOLIO RETURNS AND TARGET. 



T 



GENERATE RETURN SERIES FOR TARGET IN TERMS OF ALL ASSETS IN PORTFOLIO. | 

GENERATE CONSTRAINTS AS SECURITY TYPE, SECTOR TYPE, 
. INDIVIDUAL SECURITY WEIGHTS. 



T 



CALCULATE EACH ASSETS' COVARIANCE WITH SELECTED FINANCIAL TARGET~1 
RANK THE COVARIANCES FOR ALL POTENTIAL ASSETS OF THE PORTFOLIO. 

T 



PROCEED FROM FIRST RANKED ASSET TO END, FILLING IN WEIGHT UNTIL A 

CONSTRAINT IS REACHED, REQUESTING FOR EACH ASSET STEPWISE 
DOWN THE RANKING. 



I 



GENERATE RETURN FOR THE WEIGHTED PORTFOLIO OF ASSETS. 
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CALCULATE PARTIAL DERIVATIVE FOR EACH ASSET OF THE SELECTED 
PORTFOLIO WITH REQUEST TO THE TARGET FINANCIAL INDEX. 
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REORDER ALL ASSETS IN ORDER OF COVARIANCE AND MAKE ITERATIVE 
CHANGE IN WEIGHTS OF ASSETS TO IMPROVE OVERALL COVARIANCE. 
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REPEAT ITERATIVE iMpRuvtMtNl Oh UOVAKIANCE BY SELECTED, 
ALLOWABLE CHANGES IN ASSET WEIGHTS UNTIL CONVERGES. 
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LIABILITY RETURN 
LIABILITY BETA = 0.686 
CORRELATION = 69.77o 
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